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ABSTRACT 

We present theoretical X-ray line profiles from a range of model colliding wind 
systems. In particular, we investigate the effects of varying the stellar mass-loss rates, 
the wind speeds, and the viewing orientation. We find that a wide range of theoret- 
ical line profile shapes is possible, varying with orbital inclination and phase. At or 
near conjunction, the lines have approximately Gaussian profiles, with small widths 
(HWHM ~ O.lwoo) and definite blue- or redshifts (depending on whether the star with 
the weaker wind is in front or behind). When the system is viewed at quadrature, the 
lines are generally much broader (HWHM ~ i>oo) 5 flat-topped and unshifted. Local ab- 
sorption can have a major effect on the observed profiles - in systems with mass-loss 
rates of a few times 10~ 6 M Q yr _1 the lower energy lines (E < 1 keV) are particularly 
affected. This generally results in blueward-skewed profiles, especially when the sys- 
tem is viewed through the dense wind of the primary. The orbital variation of the line 
widths and shifts is reduced in a low inclination binary. The extreme case is a binary 
with i = 0°, for which we would expect no line profile variation. 

Key words: hydrodynamics - line: profiles - stars: binaries: general - stars: early- 
type - stars: winds, outflows - X-rays: stars 



1 INTRODUCTION 

ICherepashchukl dl97d) and IPrilutskii fc Usovl Jl97rj) pro- 
posed that the collision of the winds in a binary of early-type 
star s should produ ce copious quantities of X-ray-producing 
gas. IPollockl 1119871) found that among the 48 Wolf-Rayet 
stars observed with Einstein, the binary systems tended to 
be a few times brighter in X-rays than single stars. There 
is a simil ar trend amon g the O stars observed with Ein- 
stein JChlebowskilll989ft . This provided early observational 
evidence for colliding wind X-ray emission. 

Colliding wind binaries (CWBs) have also been studied 
extensively with numerical hydrodynamical simulations (e.g. 
i Luo. McCrav fc Mac LoyJl99d:IStevens Blondin fc Pollock! 
Il992l: iPittard fc Stevendll997t |Pittardlll998l) . These works 
investigated the broad-band X-ray properties, since the 
satellites in operation at that time had only poor spec- 
tral resolution (e.g. E/AE ~ 20 for the ASCA SIS). Nev- 
ertheless, by comparing the variable ASCA spectrum of 
7 2 Velorum with h ydrodynamical models of the system, 
IStevens et alJ lll996l) were able to confirm that the system is 
indeed a CWB, and were also able to place constraints on the 
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stars' mass-loss rates and wind velocities. This technique has 
more recently been applied to the high-resolution Chand ra 
grating spectrum of rj Carinae iPittard fc Corcoranll20o'^ . 

The unprecedented spectral resolution offered by the 
Chandra LETGS and HETGS and the XMM-Newton RGS 
(E/AE ~ 100-1000) gives us the opportunity to study in 
detail the profiles of individual X-ray lines from CWBs, 
enabling us to probe the very hot shocked gas in the 
wind- wind interaction regions. Studies of X-ray line pro- 
files from single early-type stars have already been used 
to test the standard wind-shock model of X-ray produc- 
tion, in which instabilities due to the line-driving force 
give rise to shocks distributed throughout the stellar wind 
jLucv fc Whitell98d:lLucvll982HOwocki Castor fc Rvbickl 
ll98St iFeldmeier. Puis fc Pauldracblll997ft . Line profile cal- 
culations based upon such models (with X-ray produc- 
tion effectively distributed uniformly throughout a steady- 
state or accelerating stellar wind) generally predict broad , 
blueshifted and blue-s kewed line profiles ilgnacd 1200 ll : 
lOwocki fc Cohenl l200ll) . However, blue shifted lines have 
only been observed from £ Puppis llKahn et alJ 1200 ll : 
ICassinelli et~al"ll200lD while the lines observed from S Orio- 
nis iMiller et alJl2002l) and r Scorpii dCohen et alJl2003l) are 
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far narrower than is expected from the standard wind shock 
model. 

An analogous study of X-ray line profiles from CWBs 
has not yet been made. Only relatively few observations 
of CWBs (or suspected CWBs) with high resolution spec- 
tro meters have been pu blished; the best examples are 7 2 
Vel dSkinner et alJ 1200 ill a n d n Car llCorcoran et al J 1200 it 
[ Pittard fc Corcoran! |20o3 : iLeuteneeeer. Kahn fc Ramsavl 
2003). In these works, very little attention has been paid 
to the line profile shapes. The Chandra HETGS spectrum 
of j 2 Vel is dominated by emission from H- and He-like ions 
of Ne, Mg, Si and S. The lines are broader than the expected 
thermal widths (HWHM w 500-7 00 kms" 1 ), arid the re are 
no significant wavelength shifts ISkinner et al] 12*001) . The 
lines in the Chandra HETGS spectrum of 77 Car (from H- 
like ions of Mg, Si and S and He-lik e ions of Si, S, Ca an d 
Fe) are also not significantly shifted llCorcoran et alj |2001). 
No mention is made of the widths of the lines. Similarly, 
no mention is made of line widt hs or shifts in the XMM - 
Newton RGS spectrum of r\ Car llLeuteneeeer et a""]l2005() . 
However, there is evidence for orbital variation of line shifts 
in the Chandra HETGS spectrum of WR 140 (Pollock et al., 
in prep.) - this will be discussed further in |I] 

We present here results of theoretical calculations of X- 
ray line profiles from such systems The model is described 
in g2J The calculated profiles are presented in SjU m par- 
ticular showing how the profiles vary with wind parameters 
and the viewing orientation. The results are discussed in Sj^l 
and summarized in The line profiles presented here will 
ultimately be compared with X-ray spectra of CWBs, and 
potentially offer another tool for placing constraints on the 
wind parameters of such systems, as well as providing new 
insights into the structure of the X-ray emitting regions. 



2 THE MODEL 

2.1 Hydrodynamic simulations 

The hydrodynamic simulations of CWBs were carried out 
using vh-1, an implementation of the piec ewise parabolic 
method (PPM; IColella fc Woodward! Il984h developed by 
John Blondin and co-workers at the University of Virginia 1 . 

The CWBs are characterized by the wind momentum 
ratio 1Z, given by 



1Z : 



M a va 
M b v b 



(1) 



where Ma and Mb are the mass-loss rates of the stars 
and va and vb are the terminal velocities of the winds 
iStevens et al] 1*1992'). By definition, star A (the primary) 
is the one with the more powerful wind (i.e. 1Z > 1). 

Modelling the flow in a CWB is a three-dimensional 
problem. However, except in very close binaries, the orbital 
velocity of the stars is much less than the terminal velocity 
of the winds. Since we are not modelling close binaries, Cori- 
olis forces are assumed to be negligible, and the problem re- 
duces to a two-dimensional axisymmetric flow. This greatly 
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reduces the amount of computer time needed to model the 
system. 

It is important to note that the goal of this work is to 
gain insights into where different X-ray lines originate in the 
wind-wind collision zone, and how the profiles are affected by 
different wind parameters. As a result, rather than modelling 
specific CWBs, we have chosen wind parameters that are 
representative of typical systems. In particular, the stellar 
separation is fixed at 2 x 10 13 cm (287 Rq) , mass-loss rates 
are typically ~ 10 -6 -10 -5 Mq yr _1 , and wind speeds are 
typically 2000 km s" 1 . 

For these wind parameters the winds can be assumed 
to be adiabatic, because the cooling time t coo \ is large com- 
pared with i csc , the escape time from the shock region; i.e . 
the cooling parameter \ ~ t coo i/t C sc l|Stevens et, aljfl992l) 
is large. A few simulations were carried out with radiative 
cooling included self-consistently. The inclusion of cooling 
did not significantly affect the results, and so in general it 
was not included at all. 

A further assumption is that the winds are at their ter- 
minal velocities (i><x>) when they collide. This assumption is 
valid for wide systems, since the winds of early-type stars 
typically reach their terminal velocities within a few stellar 
radii of the stellar surface. This assumption may be relaxed 
in the future for models of s pecific c l ose b inary systems, 
using the method described in iPittardl il998f) . 

Fig. is an example snapshot from one of the hy- 
drodynamical simulations, showing the temperature and 
density structure of a CWB with 1Z = v5 (Ma = 5 x 
10 -6 M yr" 1 , Mb = 1 x 10 -6 M yr" 1 , v A = v B = 
2000 km s -1 ). For these wind parameters, the cooling pa- 
rameters l|Stevens et al.lfr992T l are xa = 4.2 and xb = 11.2, 
and so the winds can indeed be assumed to be adiabatic. 
Note also that instabilities can be seen in the shocked gas, 
which may be due to either a difference in flow velocities 
along the t wo sides of the cont act discontinuity between the 
two winds dStevens et al Jl992l) or numerical effects near the 
line of centres. 



2.2 Line profile calculations 

The method used for calculating X-ray line profiles from 
the hydrodynamic simulations simply involves calculating 
the line profile for each point on the grid and summing over 
the whole grid. 

Thermal Doppler broadening is the dominant effect, 
and so each grid point is assumed to produce a Gaussian 
line profile whose height and width depend on the temper- 
ature of the gas at that point, and whose centre is blue- or 
redshifted according to the line-of-sight velocity of the gas, 



= Io exp 



-m{v - v\ os f 
2kT* 



(2) 



where P(v)dv is the intensity of the line from the ith grid 
point between a velocity shift of v and v + dv, m is the ion 
mass, T % is the gas temperature, and v\ os is the line-of-sight 
velocity of the gas. 

The height of the Gaussian 7q is calculated using the list 
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Figure 1. Density (left) and temperature (right) maps from a hydrodynamic simulation of a colliding wind binary with 1Z = \/5 
(M A = 5 x 10" 6 M Q yr -1 , M B = 1 x 10~ 6 M yr" 1 , v A = v B = 2000 km s" 1 ). 



of spectral lines available from the SPEX website 2 . This list 
gives the emissivity Q of a few thousand UV and X-ray lines 
over a range of temperatures. For an arbitrary temperature, 
the emissivity is calculated by interpolating between the tab- 
ulated values. If the temperature lies outside the range of 
temperatures for which there are data, the emissivity is as- 
sumed to be zero and the grid point is skipped over. The 
total line luminosity Li ine for the ith grid point is obtained 
by multiplying these emissivities by n\.n\ 1 V l , where n l e and 
«h are the electron and hydrogen number densities and V r 
is the volume of the emitting gas at the grid point. LJ ine is 
equal to the integral of equation @ over all v, and so 



r = g(r»^y 



2-n-kT 1 



(•3) 



Two major assumptions have to be made in order to 
use the data in the SPEX list of lines. Firstly, the plasma 
is assumed to be in collisional ionization equilibrium, and 
secondly both winds are assumed to have solar abundances. 
The latter is a valid assumption for O+O binaries, but not 
for WR+O binaries, in which the Wolf-Rayet wind has sig- 
nificantly non-solar abundances. This assumption will be re- 
laxed in the future when specific WR+O binaries are mod- 
elled. 

The line profile produced by each grid point is then at- 
tenuated by phot oabsorpt ion. The absorbing column to each 
emitting grid point is calculated by integrating the density 
numerically along the line-of-sight. The density is either ob- 
tained directly from the hydrodynamic grid, or by extrapo- 
lation when off the grid. The opacities used were calculated 
for a solar abundance plasma using xstar 3 . For a given 
line, the opacity is assumed to be constant with wavelength 
across the width of the line profile. Most of the absorption 
takes place in the cool, unshocked winds surrounding the 
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index.html 

3 http:/ /heasarc. gsfc.nasa.gov/docs/software/lheasoft/xstar/ 
xstar.html 



Table 1. The wavelengths of the Lyo lines used in this work 
(taken from SPEX list of lines). Ai and A2 are the wavelengths 
of the brighter and weaker components, respectively. Av is the 
difference in wavelengths expressed as a velocity shift (Av = 
c(A 2 - Ai)/Ai). 



ion 


Ai (A) 


A 2 (A) 


Av ( km s _1 ) 


O VIII 


18.971 


18.977 


95 


Nex 


12.132 


12.137 


124 


Mgxn 


8.419 


8.424 


178 


Si xiv 


6.181 


6.187 


291 


Sxvi 


4.729 


4.733 


254 



wind-wind collision region, which are both assumed to have 
a temperature of 3 x 10 4 K. However, opacities have also 
been calculated for higher temperatures (10 5 , 10 6 , 10 7 and 
10 8 K). The opacity of hot, shocked gas of arbitrary temper- 
ature in the wind-wind collision region is taken into account 
by interpolating between these values, though in fact this 
gas makes only a small contribution to the overall absorp- 
tion. When calculating the contribution of gas off the grid 
to the absorption, we assume its temperature is equal to the 
gas temperature at the point on the line-of-sight where we 
leave the grid. This tends to underestimate the absorption, 
as it sometimes assumes that gas off the grid is hotter than 
it actually is. The alternative is to assume that all the gas 
off the grid is cold (T = 3 x 10 4 K), which would tend to 
overestimate the absorption. However, we have found that 
there is generally very little difference in the resulting pro- 
files calculated by the two methods. 

The resulting line profile from the CWB is the sum of 
all the profiles calculated by the above method for the entire 
grid. Results for a range of CWB parameters are shown in 
the following section. 
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Figure 2. The fraction of the total line emission as a function 
of distance p from the line of centres, for the Lya lines from 
Ovin (double-dashed), Nex (dotted), Mgxil (dashed), Si XIV 
(dot-dashed) and Sxvi (solid). The mass-loss rates are Ma = 5 X 
10 -6 M Q yr _1 and M B = 1 X 10~ 6 M yr -1 . The wind speeds 
are (a) v& = vb = 2000 km s _1 , (b) va = v B = 3000 km s _1 . 
The stellar separation D = 2 X 10 13 cm. 



3 RESULTS 

Line profiles were calculated for the Lya lines from the 
abundant elements O, Ne, Mg, Si and S. These lines are 
in fact closely spaced doublets; both components were in- 
cluded in the calculations, but other nearby lines were ig- 
nored (as their emissivities are much lower than those of 
the Lya lines). The wavelengths of the lines are shown in 
Table [U 

The O to S Lya lines were chosen because they fea- 



type stars ( 


Schulz et al.ll200fi IWaldron & Cassinelli 


120011 


Kahn et alJ 


20011 ICassinelli et alJl200ll iMiller et alJ 


2002). 



They are also insensitive to the plasma density, unlike the 
fir (forbidden-intercombination-resonance) triplets from He- 
like ions. This greatly simplifies the calculations. 



3.1 Location of the line-emitting plasma 

Fig- HI shows the fraction of the total unabsorbed line emis- 
sion as a function of cylindrical distance p from the line 
of centres for each of five lines (i.e. the figure shows the 
fractional line luminosity from a series of thin cylindri- 
cal shells whose symmetry axes lie along the line of cen- 



tres). The results in Fig. |5J a ) are from a simulation with 
Ma = 5 x 10" 6 M yr -1 , M B = 1 x 10" 6 M Q yr -1 and 
va — vb = 2000 km s" 1 . The results in Fig. |2Ib) are from 
a simulation with the same mass-loss rates, but with higher 
wind velocities (3000 km s~ ). The simulations were carried 
out on an 8 x 10 13 cm square grid with a stellar separation 
D — 2 x 10 13 cm. In all cases, the values of p are normalized 
to D. 

In both cases, one can clearly see that most of the Lya 
emission from the lighter elements (e.g. O, Ne) originates 
a few times 10 13 cm from the line of centres, whereas that 
from heavier elements (e.g. S) originates much further in 
(~ 1 x 10 13 cm from the line of centres). This is exactly 
as expected, because the temperature of the plasma in the 
wind-wind collision region falls off with increasing p. In the 
inner regions it is too hot for O viii and Nex to exist in large 
amounts (these elements are mostly fully ionized), whereas 
in the outer regions it is too cool to have significant amounts 
of H-like ions of heavier elements. 

The main difference between Figs.^a) and|2Jb) is that 
for the system with higher wind velocities, the line emission 
tends to originate from larger values of p. This is because 
the higher wind velocity results in a higher post-shock tem- 
perature. 

For all the lines, the radius at which the emission peaks 
is smaller than the radius at which the plasma temperature 
equals the temperature of maximum emissivity. This is be- 
cause the emissivity Q of the lines varies fairly slowly with 
p. However, the nln^y 1 term in equation J3J varies as ap- 
proximately p~ 2 , except near the line of centres, where it 
turns over and tends to zero. As a result the emission orig- 
inates from further in than one might naively expect from 
the temperature alone. 



3.2 The profiles 

Figs. [3] to [7| show theoretical Lya line profiles for Ovin, 
Nex, Mgxil, Si xiv and Sxvi, respectively. The profiles 
have been calculated from single snapshots from the hy- 
drodynamical simulations. The model C WBs have Mb = 
1 x 10~ 6 M yr" 1 and v A = v B = 2000 km s" 1 ; the wind 
momentum ratio 1Z is adjusted by varying Ma- In all cases, 
profiles are shown for 3 values of 1Z and for 5 different view- 
ing angles 0, where 9 is the angle between the line-of-sight 
and the line-of-centres, with 9 = 0° corresponding to view- 
ing the system from behind the secondary star. For a circular 
orbit, 9 is related to the inclination i and the orbital phase 
angle 4> by 

cos 9 = sin i cos (f> (4) 

where = 0° corresponds to the secondary being in front. 
Profiles with and without absorption by the stellar winds are 
shown. The wavelength shifts are calculated relative to the 
rest wavelength of the brighter component of the Lya lines, 
and have been scaled to the wind velocity. The unabsorbed 
Ovin Lya line for 1Z = 1 has been normalized such that 
its total luminosity is equal to unity. The other lines are 
normalized relative to this line, preserving the relative line 
luminosities. 

The secondary components of the Lya lines are only 
noticeable for the shorter wavelength lines. This is because 
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Figure 3. Unabsorbcd (solid) and absorbed (dashed) O VIII Lya line profiles for a range of wind momentum ratios TZ and viewing angles 



9. The profiles are for systems with Ma = 1, 5, 9 X 10 



M 



yr 



M B = 1 x 10 _b M yr" 1 and v A = v B = 2000 km s~ 



the wavelength difference between the components (0.004- 
0.006 A; see Table Q corresponds to a larger velocity shift 
for these lines. In practice, after the lines have been folded 
with the instrumental response of the Chandra or XMM- 
Newton gratings, the two components would not be resolv- 
able. However, including both components of the Lya lines 
gives broader lines than one would get by considering only 
the brighter component. 



3.3 Variation of line parameters with viewing 
angle and wind momentum ratio 



Figs.l51to ll2l show how line parameters vary with viewing an- 
gle for each of the five lines. Each figure shows the variation 
of the mean scaled wavelength shift, the half width at half 
maximum (HWHM), and the normalized line luminosity for 
the unabsorbed and absorbed line profiles for three different 
values of 1Z (1, \fZ and 3). The horizontal dashed lines on 
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Figure 4. Unabsorbcd (solid) and absorbed (dashed) Ne X Lya line profiles for a range of wind momentum ratios 1Z and viewing angles 
6. The profiles are for systems with M\ = 1, 5, 9 X 1CT 6 Mq yr" 1 , M B = 1 X 10" 6 Mq yr~ x and v A = v B = 2000 km s _1 . 



the HWHM graphs show the thermal Doppler widths at the 
temperature of maximum line emission (for comparison). 

In the equal winds case (1Z =1), the unabsorbed lines 
are approximately symmetrical and unshifted for all view- 
ing angles. While some slight asymmetries can be seen, these 
are due to numerical effects. Theoretically, the unabsorbed 
lines for 1Z = 1 must be purely symmetrical. For 1Z > 1, 
the lines are blueshifted when the system is viewed from the 
side of the secondary (6 < 90°) and redshifted when the sys- 
tem is viewed from the side of the primary (8 > 90°). The 



maximum shifts occur when 6 = 0°, 180°. This occurs be- 
cause the shocked region is bent towards the secondary when 
7Z > 1. When the system is being viewed along the line of 
centres, the general motion of the shocked gas is either to- 
wards or away from the observer (depending on which end of 
the line of centres they are looking from), with a fairly small 
spread in line-of-sight velocities, resulting in quite narrow 
lines. However, for other viewing angles this will not be the 
case - some of the gas will be moving towards the observer, 
some will be moving tangentially to the line of sight, and 
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Figure 5. Unabsorbed (solid) and absorbed (dashed) Mgxo Lya line profiles for a range of wind momentum ratios V, and viewing 
angles 0. The profiles are for systems with Ma = 1, 5, 9 x 10~ 6 M yr -1 , M B = 1 X 10~ 6 Mq yr _1 and va = v B = 2000 km s _1 . 



some will be moving away from the observer. The result of 
this is broader lines with smaller shifts. The extreme case is 
when the system is at quadrature (0 = 90°). The unabsorbed 
lines are very broad (HWHM ~ «oo), but symmetrical and 
unshifted, apart from the presence of the secondary compo- 
nent on the redward side. This is exactly as expected, owing 
to the cylindrical symmetry of the hydrodynamic simulation 
- for every 'blob' of plasma moving towards the observer 
there's a corresponding blob on the far side of the system 
moving away at the same speed. For / 90°, the velocity 



shift (towards the blue or the red, depending on the viewing 
orientation) increases with 1Z, because the shocked region is 
bent more towards or away from the observer. 

In all cases, the HWHMs are considerably larger than 
the typical thermal widths, indicating that the line broad- 
ening is due to a wide range of line-of-sight velocities in the 
bulk motion of the hot gas, rather than merely being due to 
simple thermal Doppler broadening. 

The line luminosities increase with 1Z, because a larger 
M leads to a denser X-ray emitting plasma. However, for a 
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Figure 6. Unabsorbed (solid) and absorbed (dashed) Si XIV Lya line profiles for a range of wind momentum ratios 72 and viewing angles 
6. The profiles are for systems with M A = 1, 5, 9 x f0~ 6 M Q yr _1 , M B = 1 x 10~ 6 M Q yr _1 and v A = v B = 2000 km s -1 . 



given 1Z the unabsorbed line luminosities do not vary with 
6, since these depend on the amount of hot plasma on the 
hydrodynamical grid, which is independent of viewing angle. 

The above discussion pertains only to the unabsorbed 
line profiles. From the figures one can clearly see that ab- 
sorption has a much more significant effect on the longer 
wavelength lines, such as the O vm and Nex Lya lines. This 
is because the continuum opacity in the cold stellar winds 
increases with wavelength. The opacity at the wavelength of 
the Oviii Lya line is 27.5 times that at the wavelength of 



the Sxvi Lya line. This outweighs the fact that the longer 
wavelength lines originate from the cooler shocked gas far 
from the line of centres, and so travel through less dense 
regions of the stellar winds. 

In general, the inclusion of absorption results in 
blueward- skewed lines. This is because the gas that is mov- 
ing away from the observer is generally on the far side 
of the system, and so suffers more absorption. Indeed, for 
6 ^ 0°, 180°, the velocity shift of the absorbed line is always 
more towards the blue than the corresponding unabsorbed 
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Figure 7. Unabsorbed (solid) and absorbed (dashed) S XVI Lyct line profiles for a range of wind momentum ratios TL and viewing angle 
0. The profiles are for systems with M\ = 1, 5, 9 X 1CT 6 Mq yr" 1 , M B = 1 X 10~ 6 Mq yr~ x and v A = v B = 2000 km s _1 . 



line, though when the system is viewed from the side of the 
primary, the lines retain a mean redshift. This is because 
when the shocked region is bent away from the observer, 
there is very little blueshifted emission at all, and so the fact 
that redshifted emission is more strongly attenuated than 
blueshifted emission is less important. The range of viewing 
angles over which one expects to see redshifted emission is 
larger for the shorter wavelength lines, because absorption is 
less important for these lines. One would expect to see red- 



shifted O viii Lya emission for 9 > 140° , whereas the S xvi 
Lya line is redshifted when 6 > 90° . 

The HWHMs of the absorbed lines are generally either 
very similar to those of the corresponding unabsorbed lines, 
or considerably smaller (there is very little middle ground). 
The latter case generally occurs for the Oviii line when 
9 w 45°-75° or 6»>90°, and for the Nex line when 6»>90° 
(though the HWHMs of the absorbed and unabsorbed lines 
are again very similar as 6 approaches 180°). It also occurs 
to a lesser extent for the Mgxn line when 8 « 110°-150°. 
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Figure 8. Graphs showing how the mean wavelength shift, 
HWHM, and normalized line luminosity of the unabsorbed (solid) 
and absorbed (dashed) O VIII Lyo line profile vary with viewing 
angle. Results for 1Z = 1, \/5 and 3 are represented by crosses, 
solid squares and open circles, respectively. 
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Figure 9. Graphs showing how the mean wavelength shift, 
HWHM, and normalized line luminosity of the unabsorbed (solid) 
and absorbed (dashed) Ne X Lye? line profile vary with viewing an- 
gle. Results for 1Z = 1, \f§ and 3 are represented by crosses, solid 
squares and open circles, respectively. 



This narrowing of the absorbed lines occurs because of the 
higher opacity at the wavelengths of these lines. The profiles 
have their redward emission strongly attenuated, particu- 
larly when 1Z > 1 and the system is being viewed through 
the dense wind of the primary (9 > 90° ) . If it is attenuated to 
less than half the peak value of the blueward emission, only 
the blueshifted peak will contribute to the HWHM. This re- 
sults in a considerably smaller HWHM, though the lines do 
have large redward tails. However, it is worth noting that in 
practice such redward tails may not by observable (they may 
be lost in the underlying continuum), and so the dominant 



observational signature would be a blueshift rather than a 
broadening. On the other hand, if the redward emission is 
not that strongly attenuated, the HWHM will be approxi- 
mately equal to that of the unabsorbed line. This is generally 
the case for the Si xiv and S xvi lines. 

The luminosities of the absorbed lines decrease with 
6 for 1Z > 1. This is simply because as 9 increases, one 
changes from observing the system through the wind of the 
secondary to observing it through the denser wind of the 
primary. 
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Figure 10. Graphs showing how the mean wavelength shift, 
HWHM, and normalized line luminosity of the unabsorbed (solid) 
and absorbed (dashed) Mg XII Lya line profile vary with viewing 
angle. Results for H, = 1, V5 and 3 are represented by crosses, 
solid squares and open circles, respectively. 
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Figure 11. Graphs showing how the mean wavelength shift, 
HWHM, and normalized line luminosity of the unabsorbed (solid) 
and absorbed (dashed) Si XIV Lya line profile vary with viewing 
angle. Results for 1Z = 1, \/5 and 3 are represented by crosses, 
solid squares and open circles, respectively. 



3.4 Lines from systems with different wind 
parameters but same 1Z 

Figs. 1131 and 1141 compare how the mean scaled wavelength 
shift and HWHM vary with angle for systems with 1Z = \/5 
but different wind parameters. Results are shown for the 
O Vlll and S xvi Lya lines, which are at the two extremes of 
wavelength studied here. 

Fig. ll3l compares lines from systems with the same wind 
speeds but mass-loss rates that differ by an order of mag- 
nitude. For both lines there is no difference in the shapes 
of the unabsorbed profiles between the two systems, though 



the lines from the system with the larger mass-loss rates 
are brighter. However, the system with the higher mass-loss 
rates gives an absorbed O viii line that is narrower and more 
blueshifted. This is because this system has denser winds, 
and so the line suffers more attenuation (especially the red- 
ward emission). On the other hand the Sxvi line is not as 
strongly affected by absorption, and so there is very little 
difference between the two systems. 

For a given 1Z, the effect of increasing the mass-loss 
rates on the observed (i.e. absorbed) line luminosities is a 
combination of two competing effects. The intrinsic line lu- 
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Figure 12. Graphs showing how the mean wavelength shift, 
HWHM, and normalized line luminosity of the unabsorbed (solid) 
and absorbed (dashed) S XVI Lyct line profile vary with viewing 
angle. Results for H, = 1, V5 and 3 are represented by crosses, 
solid squares and open circles, respectively. 



minosity increases as M , whereas the column density in- 
creases as M. This means that the observed line luminos- 
ity scales as M 2 e~ aM , where a is a constant that depends 
on 8 and the opacity for the line in question. For the sys- 
tems described here, we find that the observed line lumi- 
nosity always increases with increasing M (for all lines and 
for all 8). However, for larger mass- loss rates (few times 
1(T 5 M yr" 1 ) the increasing absorption starts to domi- 
nate, and when 8 ~ 180° (which is when the column density 
is largest) the observed luminosities of the longer wavelength 
lines start to decrease with increasing M. 



Fig. 1141 compares the same two lines from systems with 
the same mass-loss rates but different wind velocities. Once 
the wavelength shift and the HWHM have been scaled to the 
wind velocity, there is very little difference between the two 
systems. The O Vlll line is generally slightly more blueshifted 
for the system with the lower wind velocity (since, for the 
same mass-loss rate, this results in a denser wind and so 
more absorption). On the other hand, the wavelength shift 
of the S xvi line (whether to the blue to the red) is smaller 
for the system with the lower wind velocity. This may be 
because the lower wind speed results in a lower post-shock 
temperature, and so the line forms nearer the line of cen- 
tres for this system. The shocked region is less bent towards 
(away from) the observer near the line of centres, resulting 
in a smaller blueshift (redshift) as 8 approaches 0° (180°). 



3.5 Time variability 

Eddies/swirls in the shocked gas are frequently seen in the 
hydrodynamical simulations of the CWBs used for this work. 
These travel from the line of centres to the edge of the grid 
(typically 3.4 x 10 13 cm from the line of centres) on a time- 
scale of ~ 200 ks. These swirls may be due to either in- 
stabilities alon g the conta ct discontinuity between the two 
winds llStevens et alll992T) or numerical effects near the line 
of centres. 

In order to investigate whether these swirls have any 
significant effect on the results, we have compared line pro- 
files generated from snapshots of the hydrodynamical simu- 
lations separated by ~ 25 ks. Over the course of 10 such 
snapshots, we find that the calculated wavelength shifts, 
HWHMs and line luminosities typically vary by only a few 
per cent. This suggests that any observed time variability 
in the line profiles will be due to changes in the viewing 
orientation (because of the orbital motion), rather than to 
instabilities, etc. in the shocked gas. 



3.6 Simulations of observed profiles 

We have carried out some Monte Carlo simulations of ob- 
served line profiles from CWBs. Firstly, a theoretical pro- 
file is convolved with a Gaussian, which is a good first or- 
der approximation to the Chandra HETGS line response 
function (Chandra Proposers' Observatory Guide 4 , ver. 5.0, 
§8.2.2). The width of the Gaussian used is appropriate for 
the wavelength of the line in question 5 . This convolved pro- 
file is then used as a probability distribution function in 
the Monte Carlo simulations. For the purposes of these sim- 
ple calculations, the background/continuum emission in the 
vicinity of the line is assumed to negligible. 

Fig. 1151 shows the results of simulated observed Mg XII 
Lya line profiles from a CWB with 1Z — yE for three dif- 
ferent viewing angles (8 = 0°, 45°, 90°). In each case there 
are assumed to be 200 photons in the line. The results have 
been binned up to 0.01 A. These results indicate that the 
features described in H3.3l should easily be detectable by the 



4 http:/ /asc. harvard.edu/proposer/POG/html/ 

5 The line resolving power data are taken from the 
file hetghegD1996-ll-01res_conN0002.rdb, available from 
http: / / space.mit.edu/HETG / caldb / ard_stat_hetg.html 
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Figure 13. Graphs showing how the mean wavelength shift and HWHM of the unabsorbed (solid) and absorbed (dashed) Ovin and 
S XVI Lya lines vary as a function of viewing angle for two different systems with TZ = V5. The solid squares show results for a system 
with Ma = 5 X 10 -6 Mq yr -1 , Mb = 1 x 10~ 6 Mq yr -1 and v\ = v b = 2000 km s _1 . The crosses show results for a system with 
M A = 5 X 10 -7 M Q yr -1 , M B = 1 X 10 -7 Mq yr" 1 and v A = v B = 2000 km s _1 . 



Chandra HETGS: near conjunction (9 = 0° or 180°) the 
lines are narrow with measurable blue- or redshifts, whereas 
near quadrature (0 = 90°) the profiles are very broad and 
fairly fiat-topped with negligible wavelength shifts. Note 
that for this system va = vb = 2000 km s _1 , corresponding 
to AA = 0.056 A. Hence, near quadrature the HWHM of 
the observed line is an appreciable fraction of the wind ve- 
locity, as expected. Note also that the photon statistics and 
the instrumental response can result in an observed profile 
that is distinctly different from the theoretical profile for 
low counts: compare the middle panel of Fig. 1151 with the 
appropriate profile in Fig.|S] 



4 DISCUSSION 

We have seen that X-ray lines from CWBs show significant 
variation of wavelength shift and width as a function of view- 
ing orientation, which is itself a function of orbital phase and 
inclination. The observation of phase-locked variation in the 
X-ray line profiles would therefore provide convincing evi- 
dence for X-ray emission from a colliding wind region. In 
this respect, an eclipsing binary would be the most inter- 
esting system, since during the course of its orbit the full 
range of viewing angles discussed above would be observed. 
In contrast, a system with an inclination of 0° would always 
be observed at 9 = 90°. Any phase-locked variation in this 



case would just be associated with a change in orbital sepa- 
ration and hence the density of the shocked gas. This would 
probably just manifest itself as a change in line luminosity, 
not as a change in line shift or width. 

Line profiles that vary during the orbit of a CWB have 
already been observed in the Chandra spectra of WR 140 
(Pollock et al., in prep.). The first observation was taken 
just before periastron at an orbital phase o f 1.982 accord- 
ing to the ephemeris of lWilliams et all dl990l) . At this phase 
the wind coll ision region is bent appr oximately towards the 
observer (see lWhite fc Becked lll995t) for a diagram of WR 
140's orbit), and the lines are significantly blueshifted, as 
would be expected. The second observation was taken just 
after periastron at a phase of 2.027. Although there are con- 
siderably fewer counts in this second spectrum, the lines are 
clearly much broader, and there is also some evidence of red- 
shifts. While we have not compared in detail our theoretical 
profiles with the spectra of WR 140, these observations seem 
to provide at least a qualitative confirmation of some of our 
general predictions. 

A more detailed comparison of our model line profiles 
with observed high-resolution X-ray spectra of CWBs may 
enable us to place constraints on certain parameters of the 
system. In practice, one would just fit to the shifts and 
widths of the lines, and scale the luminosities of the the- 
oretical lines to the observed luminosities. (This approach 
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Figure 14. Graphs showing how the mean wavelength shift and HWHM of the unabsorbed (solid) and absorbed (dashed) O VIII and 
S XVI Lya lines vary as a function of viewing angle for two different systems with TZ = V5. The solid squares show results for a system 
with Ma = 5 X 10 — 6 Mq yr _1 , Mb = 1 x 10 -6 Mq yr -1 and da = v b = 2000 km s — 1 . The open circles show results for a system with 
lit A = 5 x 10~ 6 
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of fitting the line shapes, not the luminosities , has already 
been used bv lKramer. Cohen fc Owockil ll2003t) to model the 
line emission from the single O star £ Pup.) For a CWB, 
the main parameters would be the wind momentum ratio 
1Z and the viewing angle 9, although it may also be possi- 
ble to place constraints on the individual wind parameters. 
For this method to be effective, one would want to fit as 
many lines as possible simultaneously. While we have only 
presented results for Lya lines here, the method described 
could be used for any line that is not affected by the density 
(e.g. Ly/3 and Ly7 lines, and the resonant component of the 
fir triplets from He-like ions). 



The model presented here for calculating X-ray line pro- 
files from CWBs is fairly simple, but it contains the essen- 
tial physics and gives an overview of how the profiles vary 
with wind parameters and viewing angle. When the model 
comes to be applied to real systems, extra physics will be in- 
cluded in the calculations as is appropriate for the systems in 
question. This could include modifying the hydrodynamical 
simulations (e.g. radiative cooling for systems with small or- 
bital separations and/or large mass- loss ratesl|Steven^3^iL| 
1992) or radiative driving for close systems iPittar di ll 998)). 
or modifying the actual line profile calculations (e.g. taking 
into account non-solar abundances for WR+O binaries). 



5 SUMMARY 

We have presented theoretical X-ray line profiles from model 
colliding wind systems for a range of mass-loss rates, wind 
velocities and viewing angles. A wide range of profile shapes 
is possible, varying with orbital inclination and phase. 

For a high inclination binary with 1Z > 1 there is a 
general blue/redshift at phases and 0.5 (i.e. at the two 
conjunctions). The profiles are narrow (HWHM ~ O.luoc) 
and approximately Gaussian. By contrast, at quadrature 
the profiles are generally very broad (HWHM ~ «oo), flat- 
topped and unshifted. Local absorption has a major effect 
on the observed profiles. When a system with relatively large 
mass-loss rates (few times 10~ 6 Mq yr -1 ) is viewed through 
the dense wind of the primary, the lines are generally skewed 
towards the blue. This is because the redshifted emission is 
strongly attenuated. This effect is most noticeable for the 
lower energy lines (i5 < 1 keV). These blueward-skewed pro- 
files generally have large redward tails. In practice these tails 
may not be observable, and so the line would be seen as fairly 
narrow and blueshifted, rather than broad and blueward- 
skewed. 

A binary with a low orbital inclination will exhibit less 
orbital variation of line width and shift, the extreme case 
being a binary with i = 0°, for which we would expect no 
orbital variation in the line profile. 

After the line shifts and widths have been scaled to 
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Figure 15. Simulated Mgxn Lya profiles as would be observed 
by the Chandra HETGS-HEG. The profiles are for a system with 
It = VI (M A = 5 X 10" 6 M Q yr" 1 , M B = 1 X 10~ 6 M yr"\ 
= = 2000 km s -1 ) observed from three different viewing 
angles (9 = 0°, 45°, 90°). There are 200 photons in the line in 
each case. The dashed vertical line shows the rest wavelength of 
the brighter component of the line (see Table HI . 



the wind velocity, there is very little difference between the 
results from systems with va ~ vb ~ 2000 km s _1 and with 
va = vb = 3000 km s - . The former systems have denser 
winds for the same mass-loss rate, and hence absorption has 
a larger effect. However, it does not significantly affect the 
results. 

Simple Monte Carlo simulations of observed X-ray line 
profiles show that the features discussed above will be de- 
tectable in good Chandra HETGS spectra of CWBs. Com- 



paring our theoretical profiles with such spectra potentially 
offers another method for determining the wind parameters 
of CWBs, and may also give us new insights into the struc- 
ture of the wind- wind collision regions. 
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